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Chapter 1 

/z-transforms and orthogonal polynomials 


Dominique Bakry and Olfa Zribi 


1.1 Introduction 


When the first author of this paper was a student, he was attending the DEA course of Marc Yor, about Brownian 
motions and the many laws that one would compute explicitly for various transformations on the trajectories. It 
looked like magic, and was indeed. In particular, the fact that conditioning a real Brownian motion to remain 
positive would turn it into a Bessel process in dimension 3, that is the norm of a 3-dimensional Brownian motion, 
seemed miraculous. Of course, there are much more striking identities concerning the laws of Brownian motion 
that one may find in the numerous papers or books of Marc Yor (see [26] for a large collection of such examples). 
The same kind of conditioning appears in many similar situations, and specially in geometric models. This is 
related to the fact that we then have explicit h (or Doob)- transforms. 

This relation between conditioning and /z-transform was first put forward by J.L. Doob [11], and is described in 
full generality in Doob’s book [12]. However, this kind of conditioning has been extended in various contexts, and 
very reader friendly explained by Marc Yor and his co-authors, in particular in [29, 27]. The fact that conditioning 
a d-dimensional model to remain in some set produces a new model in the same family (whatever the meaning of 
"family"), moreover with dimension d + 2, appears to be a general feature worth to be further understood. It turns 
out that the most known models have a common explanation, due to an underlying structure related to orthogonal 
polynomials. The scope of this short note is to shed light on these connections. 

The paper is organized as follows. In Section 1.2, we present the langage of symmetric diffusion operators that 
we shall use in the core of the text, and explain what /z-transforms are. Section 1.3 gives a few classical and known 
examples (some of them less well known indeed). They all follow the same scheme, explained in Section 1.4, 
which provides the general framework, related to the study of orthogonal polynomials which are eigenvectors of 
diffusion operators. The last Section 1.5 provides further examples, as applications of the main result, inspired 
from random matrix theory. 
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1.2 Symmetric diffusion operators, images and and /i-transforms 
1.2.1 Symmetric diffusion operators 

We give here a brief account of the tools and notations that we shall be using throughout this paper, most of them 
following the general setting described in [2]. A symmetric diffusion process (f,) on a measurable space E may 
be described by its generator 2zf, acting on a good algebra ,<// of real valued functions (we shall be more precise 
about this below). The diffusion property is described through the so-called change of variable formula. Namely, 
whenever f = (ft,-- - , G &/ p , and if <£ : R p 1-4 R is a smooth function such that <P(f) G srf together with 
<?,<£(/) and <?,;<£(/), ViJ = 1 • • • n. then 

^W)) + (1-D 

i ij 

where F(/,g) is the square field operator (or carre du champ), defined on the algebra .cF through 

r(f,g) = 

This change of variable formula (1.1) is some "abstract" way of describing a second order differential operator 
with no 0-order term. It turns out that the operators associated with diffusion processes satisfy F(/,/) > 0 for any 
/ G srf, and that the operator F is a first order differential operator in each of its argument, that is, with the same 
conditions as before, 

rw),g)) = (i.2) 

i 

In most cases, our set E is an open subset LI C R", and the algebra ,cF is the set of smooth (that is functions 
n i-4 R. Then, using formula (1.1) for a smooth function / : Cl h 4 R instead of <J> and (xi,--- ,x„) instead of 
(/i, • • • ,/„), we see that Jz? may be written as 

*{f) = I g ii {x)df j f + Y d b\x)d i f, (1.3) 

ij i 


and similarly 

r(f,g) = Y J s ll {x)d i fd j g. 
ij 

In this system of coordinates, = F(x/,xy) and b‘ = J5f(x ; ). The positivity of the operator F just says that the 
symmetric matrix (g' J )(x) is non negative for any x G LI , which is usually translated into the fact that the operator 
is semi-elliptic. In the same way, the absence of constant term translates into the fact that for the constant function 
1, that we always assume to belong to the set srf, one has 2z?( 1) = 0, which is an easy consequence of (1.3). 

It is not always wise to restrict to diffusion operators defined on some open subsets of R". We may have to 
deal with operators defined on manifolds, in which case one may describe the same objects in a local system of 
coordinates. However, using such local system of coordinates in not a good idea. In Section 1.5.1 for example, we 
shall consider the group SO(d) of tZ-dimensional orthogonal matrices. The natural algebra ,e/ that we want to use 
is then the algebra of polynomial functions in the entries (/«, ,) of the matrix, and the natural functions <f> acting on 
it are the polynomial functions. Since the polynomial structure will play an important role in our computations, it 
is not wise in this context to consider local system of coordinates (the entries of the matrix cannot play this role 
since they are related through algebraic relations). 

Coming back to the general situation, the link between the process (f,) and the operator 2z? is that, for any 

/G^, 



/(&)-/(§>)- f &{f)&)ds 
Jo 

is a local martingale, and this is enough to describe the law of the process starting from some initial point tgo=xG 
E , provided the set of functions ■<// is large enough, for example when ■<// contains a core of the so-called domain 
of the operator Jzf, see [2], chapter 3, for more details. 

The law of a single variable <* ( , when go — x, is then described by a Markov operator P t , as 

Pt(f)(x)=E x (f^t)), 

and, at least at a formal level, P t = explt.bf) is the semigroup generated by ££. 

In most of the cases that we are interested in, the operator Jz? will be symmetric in some L 2 (/i) space. That is, 
for some subset sz/q of ///, which is rich enough to describe P, from the knowledge of 2z? (technically, as mentioned 
above, a core in the domain $>{££)), one has, for f,g in ,c^o 


This translates into the integration by parts formula 

J fj?( g )dp = -Jr(f, g )dn. (i.4) 

For an operator given in an open set £2 C M" by the formula (1.3), and when the coefficients g ,J and b' are 
smooth, one may identify the density p (x) of the measure ii, when p (x) > 0, by the formula 

-W) = E WA/). 

which gives 

b' = Y^(g ,J djlogp+djg ,J ), (1.5) 

j 

an easy way to recover p up to a multiplicative constant provided (g'i) is non degenerate, that is when 2zf is elliptic. 
We call this measure p the reversible measure. Indeed, whenever the measure p is a probability measure, and under 
this symmetry property, then the associated process (q ) has the property that, whenever the law of go is p, then 
for any t > 0 the law of G [0, f]) is identical to the law of (2g s ,s € [0, f]). This justifies in this case the name 

"reversible", which we keep in the infinite mass case, following [2], 

Through the integration by parts formula, the operator 2z? (and therefore the process and the semigroup them¬ 
selves, provided we know something about a core in the domain), is entirely described by the triple (fl.r.p), 
called a Markov triple in [2], 

Thanks to the change of variable formula (1.1), it is enough to describe an operator in a given system of coor¬ 
dinates (x‘) to describe 2z?(x l ) = b' and r(x',x J ) = g l K Indeed, this determines 2z?(<J>(x ! )), for any <f> at least 
As outlined earlier, we do not even require that these functions x‘ form a coordinate system. They may be redun¬ 
dant (that is more variables than really necessary, as for example in the SO(d) mentioned above), or not sufficient, 
provided the computed expressions depend only on those variables, as we do for example in Section 1.5. 

Moreover, it may be convenient in even dimension to use complex variables, that is, for a pair (x,y) of functions 
in the domain, to set z =x + iy and describe Jz?(z) = Jz?(x) + iJz?(y), r(z,z ) = E(x.x) — T(y,y) +2ir(x,y) and 
r(z,Z ) = r (x.x) ■ r(y.y), and similarly for many pairs of real variables, or a pair of a real variable and a complex 
one. This will be used for example in paragraphs 1.3.4 and 1.5.2. However, we shall be careful in this case to apply 

only to polynomial functions in the variables (x,y), replacing x by j(z + z) and y by y (z —z). Then, the various 
change of variable formulae (on «£? and F) apply when considering z and z. as independent variables. 
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As we already mentioned, it may happen that we can find some functions A ,.i l,--- ,k such that, for any i, 
66(X,) depend only on (Ai,--- ,Xf) and that the same is true for r(Xj,Xj) for any pair (i,j). Then, writing X = 
(X ], • • • ,Xk) G R k , setting B' (A) = 66 (X 1 ) and G‘ J (A ) = T (A,-, A/), one writes for any smooth function <P : R* >->• R, 
66(<P(X)) = 6f(&)( A), where 

^ = Y l G ij (X)d? j + ^B i (X)d i , 
ij i 

which is a direct consequence of formula (1.1). When such happens, the image of the process (| f ) with generator 
, 16 ' under the map A is again a diffusion process (q ) with generator , 16 '. In this situation, we say that , 16 ’ is the 
image of 66 through the map A. 

Some caution should be taken in this assertion concerning the domains of the operators, but in the examples 
below all this will be quite clear (our operators will mostly act on polynomials). When 66 is symmetric with respect 
to some probability measure li, then 66 is symmetric with respect to the image measure il of if through A. With 
the help of formula (1.5), it may be an efficient way to compute ii. 


1.2.2 h-tranforms 

Given some diffusion operator 66 on some open set in R^, we may sometimes find an explicit function h, defined 
on some subset 61\ of 61, with values in (0,°°) such that 66(h) = Xh, for some real parameter A > 0. We then look 
at the new operator 66 > ( lt \ acting on functions defined on Q\, described as 

is another diffusion operator with the same square field operator than 66 . This is the so-called h (or Doob’s) 
transform, see [11, 12, 2], Indeed, thanks to formula (1.1), one has 

66 (h \f) = 66(f) + 2 r{\ogh,f). 

When the operator 66 is symmetric with respect to some measure li, then 66 ih) is symmetric with respect to 
dfJ,i t = h 2 djl. 

Considering functions with support in 6l\, the application D is an isometry between L 2 (/i/ ; ) and L 2 (/i). 

It is worth to observe that 66^ = D 1 (66 — Ald)D : every spectral property (discreteness of the spectrum, cores, 
etc.) is preserved through this transformation. 

For example, if / G L 2 (/i) is an eigenvector of 66 with eigenvalue —A], then f/h is an eigenvector of 66^ h> with 
eigenvalue — (Ai +A). 

Also, at least formally, for the semigroup p}^ associated with 2z?W, one has 

P W(f) =e - Xt jP t (hf). 

In general, one looks for positive functions h which vanish at the boundary of 61\, and there is a unique such 
function h satisfying 66(h) = —Xh, usually called the ground state for 66 on 61\. This situation appears in general 
in the following context. When 66 is elliptic on 61 C R", and whenever 61\ is bounded, with fl\ C 11, there one 
may consider the restriction of 66 on 61\. If we impose Dirichlet boundary conditions, then the spectrum of this 
operator consists of a discrete sequence 0 > Ao > Ai > • • • > A„ • • •. The eigenvector h associated with Ao is strictly 
positive in 61\ and vanishes on the boundary d61\. This is the required ground state h of the operator on 61 \. 



In probabilistic terms, the operator W is the generator of the process (q,), conditioned to stay forever in the 
subset 12 1 . However, this interpretation is not that easy to check in the general diffusion case. We shall not be 
concerned here with this probabilist aspect of this transformation, which is quite well documented in the literature 
(see [12] for a complete account on the subject, and also [29, 27] for many examples on conditioning), but rather 
play around some algebraic aspects of it in concrete examples. However, for the sake of completeness, we shall 
briefly explain the flavor of this conditioning argument in the simplest example of finite discrete Markov chains, 
where all the analysis for justification of the arguments involved may be removed. 

For this, let us consider a finite Markov chain (X„) on some finite space E, with probability transition matrix 
P{x,y ), (x,y) £ E 2 , which would play the role of P\ in the diffusion context. For simplicity, let us assume that 
P{x,y) > 0 for any (x,y) £ E 2 . Consider now a subset A £ E, and look at the restriction P\ of the matrix P to A x A. 
The Perron-Frobenius theorem asserts that there exists a unique eigenvector Vo for Pa, associated with a positive 
eigenvalue /i 0 , which is everywhere positive. This eigenvector Vo corresponds to the ground state h described above 
in the diffusion context. Then, one may look at the matrix Q on ,4 x A, defined through 


Q{x,y) 


Vo(y) 

MoVo(x) 


P(x,y), 


which is a Markov matrix on A x A. This Markov matrix Q plays on A the role of exp(22f ( ^) when h is the ground 
state on Q. \. 

Fix now n > 0 and N > n. Let A ,y be the event (A)] £ A,-- - , Xa £ A). For the initial Markov chain (X„) 
with transition matrix P and for Xq = x £ A, consider now the law of (Xq, ■ ■ ■ ,X n ) conditioned on Aj y. When 
F(Xq ,• • • ,X n ) = fo(X q)■■■ f n (X n ), it is quite easy to check that 


E(F(Xo,---X„) 1 An ) 
E(1 a n ) 


1 i? 

Q N (l/V 0 )(x) E 


(F{X 0 ,---,X n )Q N - n {\/Vo){X n )), 


where E denotes the expectation for the law of a Markov chain with matrix transition Q. 

Now, using the irreducibility of the Markov matrix Q, one sees that, when N goes to infinity, both Q N ~ n ( \/Vo){X n ) 
and Q n ( 1 /Vq)(x) converge to f pdv, where v is the (unique) invariant measure for the matrix Q. In the limit, we 
recover the interpretation of the transition matrix transition Q as a matrix of the conditioning of the Markov chain 
(X„) to stay forever in A. 

Coming back to the general case, it is worth to observe that, at least formally, the transformation JSf i —y 2S?W is 
an involution. Indeed, Sf^(l) = —j and (_£f( ,2 ))( 1 / ft ) = jzf. However, in the usual context of ground states, the 
interpretation of the associated process as a conditioning is more delicate, since 1 /h converges to infinity at the 
boundary of the domain fl\. 

It is not in general easy to exhibit explicit examples of such ground states h, but there are many very well known 
examples in the literature. We shall show that in the realm of diffusion processes which are associated to families 
of orthogonal polynomials, there is a generic argument to provide them, and that this family of examples cover 
most of the known ones, either directly, either as limiting cases. 


Remark 1. Observe that, beyond the case where h is a positive eigenvector for «£?, one may use the same transfor- 
mationforany positive function/;. One may then look at (/) = = 2z?(/) +2F (log/;,/) +Vh, where 

V = In particular, with operators in R" of the form St ?(/) = A (/) + Vlog V • V/, which have reversible 

measure Vdx, one may use h = V '^ 2 , which transforms in an isospectral way AX' into a Shrodinger type operator 
Af + Vf, associated with Feynman-Kac semigroups. This allows to remove a gradient vector field, the price to 
pay is that one adds a potential term. This technique may be used to analyse spectral properties of such symmetric 
diffusion operators through the techniques used for Shrodinger operator (see [2], for example). 
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1.3 Some examples 
1.3.1 Bessel operators 

We start form the Brownian motion in R. The operator Jz? is given by Jz?(/) = \ f". Here, F(/,/) = j /' 2 and /r 
is the Lebesgue measure. If we consider Q = (0,°°) and h = x, one has A = 0 and Jf' W(/") = j(/" + |/'). This 
last operator is a Bessel operator SB 3 . More generally, a Bessel process Bes{n) with parameter zz has a generator in 
(0,°°) given by 

2 x 

and it is easily seen, when zz > 1 is an integer, to be the generator of \\B t ||, where ( B,) is an zz-dimensional Brownian 
motion (indeed, S$ n is the image of the Laplace operator \A under x ||x||, in the sense described in Section 1.2). 
This SB 3 operator is also the generator of a real Brownian motion conditioned to remain positive. Observe however 
that the function h is this case does not vanish at the infinite boundary of the set (0,°°), and that the probabilistic 
interpretation would require some further analysis than the one sketched in the previous section. 

From formula (1.5), it is quite clear that a reversible measure for the operator SB n is x "~ 1 dx on (0,°°), which for 
n £ N*, is also, up to a constant, the image of the Lebesgue measure in R" through the map x *-> ||x||. 

This /z-transform may be extended to the general Bessel operator. Indeed, for any zz > 0, one may consider the 
function h„(x) = x 2 ~”, for which SB n {h n ) = 0, and then SBn '"' 1 = SB^,,. 

The change of SB n into SB^ n is perhaps more clear if we consider the generator through the change of variable 
x 1 —^ x 2 , that is if we consider the generator of the process (<i; 2 ) instead of the process (<y) with generator SB n . A 
simple change of variable provides the image operator 

& n (f)= 2 xf + nf, (1.6) 

for which the reversible measure has density p (x) = x^ n ~ 2 ^ 2 , and the function h is nothing else than 1 /p. 

Under this form, we shall see that is a particular case of a phenomenon related to orthogonal polynomials, 
developed in Section 1.4, although here there are no polynomials involved here, the reversible measure being 
infinite. 

Remark 2. It is not hard to observe that for 0 < n < 2, the process (£ f ) with associated generator SB n , and starting 
from x > 0 reaches 0 in finite time. Then, S& 4 -„ is the generator of this process conditioned to never reach 0. 
However, it is well known that the Bessel operator is essentially self-adjoint on (0,°°) as soon as n > 3 (see [2], 
page 98, for example). This means that the set of smooth function compactly supported in (0,°°) is dense in the 
L 2 domain of SB n . Since this is a spectral property, it is preserved through /z-transform and this also shows that it 
is also essentially self adjoint for any n < 1. In particular, there is a unique symmetric semi-group for which the 
generator coincides with SB n on the set of smooth compactly supported functions. On the other hand, for 1 < zz < 2, 
since the associated operator hits the boundary in finite time, there are at least two such semigroups with as 
generator acting on smooth functions, compactly supported in (0,°°) : the one corresponding to the Dirichlet 
boundary condition, corresponding to the process killed at the boundary {x = 0}, and the one corresponding to 
the Neuman boundary condition, corresponding to the process reflected at the boundary. Through /z-transforms, 
one sees then that there are also at least two positivity preserving semi groups in the case 2 < zz < 3, which may 
be a bit surprising since then the associated process does not touch the boundary. However, although the Dirichlet 
semigroup is Markov (P f (l) < 1 ), its /z-transform is Markov (/) ( 1 ) = 1 ), while the /z-transform of the Neuman 
semigroup (which is Markov), satisfies P f ( 1 ) > 1 . 



1.3.2 Jacobi operators 


This is perhaps the most celebrated case of known explicit /z-transform, since it is closely related in some special 
case to the Fourier transform on an interval. The Jacobi operator on the interval (—1,1) has generator 

A, p if) = ( J - ■* 2 )/" - ((«+P )* -+« - P )/' 

and is symmetric with respect to the Beta measure on (—1,1) which is C a p(l — x) a ~ l (l + x)^~ l dx, C a p being 
the normalizing constant. We always assume that a. p >0. There is a duality through /z-transforms exchanging 
p and J? 2 -a 2 -/ 3 * the function h being (1 —x) 1_ “(l — x) 1- ^, that is, as in the Bessel case in the appropriate 
coordinate system, the inverse of the density measure. 

In a similar way that the Bessel process may be described as a norm of a Brownian motion, one may see the 
symmetric Jacobi operator (a = p) as an image of a spherical Brownian motion in dimension 2a. Namely, if one 
considers the unit sphere §” in E" 1 1 , and looks at the Brownian motion on it (with generator A§n being the Laplace 
operator on the sphere), and then one looks at its first component, one gets a process on (—1,1) with generator 
Jzf«/2,n/2_ (We refer to paragraph 1.3.5 for details about the spherical Laplacian, from which this remark follows 
easily, see also [2, 31]). One may also provide a similar description in the asymmetric case, when the parameters a 
and P are half integers. In this case, 2£ a p is, up to a factor 4, the image of the spherical Laplace operator acting on 
the unit sphere § 2 “+ 2 P- 1 through the function X : § 2 a + 2 P - 1 ^ [—1,1] defined, forx = (xj, • • • ,X 2 a + 2 p) G R 2a+2 ^ 
as 

2a 

X(x) = -l+2 £x 2 . 

;=i 

The operator ^ a p may be diagonalized in a basis of orthogonal polynomials, namely the Jacobi polynomials. 
They are deeply related to the analysis on the Euclidean case in the geometric cases described above. For example, 
when a = P is an half-integer, then, for each degree k, and up to a multiplicative constant, there exists a unique 
function on the sphere which depends only on the first coordinate and which is the restriction to the sphere of 
an homogeneous degree k harmonic polynomial in the corresponding Euclidean space : this is the corresponding 
degree k Jacobi polynomial (see [31, 2] for more details). In other words, if I\ (x) is one of these Jacobi polyno¬ 
mials with degree k corresponding to the case a = p = n/2, then the function (xi, • • • ,x„+i) >->• ||x|| A(imr) is an 
homogeneous harmonic polynomial in R ,,+ 1 . A similar interpretation is valid in the asymmetric case, whenever the 
parameters a and p are half-integers, if one reminds that the eigenvectors of the Laplace operator on the sphere 
are restriction to the sphere of harmonic homogeneous polynomials in the ambient Euclidean space (see [31]). 

For a = P = 1 /2, r J? a p this is just the image of the usual operator f" on (0, n) through the change of variables 
0 i —y cos(O) = x. More generally, in the variable 0, may be written as 

a d 2 (a + P — l)cos(0) + a — P d 
A,/3 = ^2 + sin(0) c/0 

For a = p = 1/2, corresponding to the arcsine law, the associated orthogonal polynomials p)/ 2A ~ are the 
Chebyshev polynomials of the first kind, satisfying 

p„ 1/2 ’ 1/2 (cos(0)) =cos(zz0). 

For a = P =3/2, corresponding to the semicircle law, they correspond to the Chebyshev polynomials of the 
second kind, satisfying the formula 


sin(0)/ ,2/, “’ 3//2 (cos(0)) = sin (zz0). 
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These formulae indeed reflect the /z-transform between ^V 21 / 2 and ^ 3 / 2 > 3 / 2 . While /j] 22 ' /2 (cos(0)) is a basis 
of L 2 ((0, 7t),dx) with Neuman boundary conditions, sin(0)f’,f /,2;3//2 (cos(0)) is another basis of L 2 ((0, n),dx), 
corresponding to the Dirichlet boundary condition. This is the image of the eigenvector basis for j/f 3 / 2 - 3 / 2 through 
the inverse h transform, the function h being in this system of coordinates nothing else than (sin0) . 

For n = 1, one gets the projection of the Brownian motion on the circle, which is locally a Brownian motion 
on the real line, up to a change of variables. The first coordinate x\ on the sphere plays the role of a distance to 
the point (1,0,• • • ,0) (more precisely, arccos(xi) is the Riemannian distance on the sphere from (1,0,• ,0) to 
any point with first coordinate x{), and we have a complete analogue of the case of the one dimensional Brownian 
motion. Namely, 

Proposition 1. The Brownian motion on the half interval (identified with the circle) conditioned to never reach the 
boundaries is, up to a change of variable, the radial part of a Brownian motion on a 3 dimensional sphere. 


1.3.3 Laguerre operators 

This is the family of operator on (0,°°) with generator 

-£(«) (/) = xf + (a - x)f, 
which is symmetric with respect to the gamma measure 

dp 1 ' 01 ' 1 = C a x a ~ 1 e~ x dx. 

For a > 0, the Laguerre family of operators is another instance of diffusion operators on the real line which may 
be diagonalized in a basis of orthogonal polynomials : these polynomials are the Laguerre polynomials, and are 
one of the three families, together with Jacobi polynomials and Hermite polynomials, of orthogonal polynomials 
in dimension 1 which are at the same time eigenvectors of a diffusion operator, see [3], The Laguerre operator is 
closely related to the Ornstein-Uhlenbeck operator defined in (1.7), and plays for this operator the same role that 
the one played by Bessel operators for the Euclidean Brownian motion. 

It is indeed quite close to the Bessel generator under the form (1.6), and in fact the Bessel operator may be 
seen as a limit of Laguerre operators under proper rescaling. It is also a limit of asymmetric Jacobi operators, also 
under proper rescaling (see [2]). The function h = x l ~ a satisfies 2z?( a )(/z) = (a — l)/z, and the /z-tranform of 

is ^(2-a)- 

As mentioned above, when a is a half-integer n /2, the Laguerre operator may be seen as the radial part of the 
Ornstein-Uhlenbeck operator in K" with generator 

f£ ou = A —xV, (1.7) 

which is symmetric with respect to the standard Gaussian measure. More precisely, for a = n/2, = 

2^«5f( a )/j It is therefore an image of the n -dimensional Ornstein-Uhlenbeck operator in the sense of Sec¬ 

tion 1.2. In other words, the Laguerre process with generator 2«5f(„y 2 ) is nothing else than the squared norm of an 
Ornstein-Uhlenbeck process in K". For a = 1/2, this corresponds to the modulus of a one dimensional Ornstein- 
Uhlenbeck, that is the one dimensional Ornstein-Uhlenbeck operator itself on (0,°°), and we get, as the particular 
case for n = 1 /2, 

Proposition 2. The law of an Ornstein-Uhlenbeck operator in dimension 1, conditioned to remain positive is the 
same as the law of the norm of a 3-dimensional Orntein-Uhlenbeck operator. 



1.3.4 An example in M 2 

The following example, less well known, had been pointed out by T. Koornwinder [22], not exactly under this form 
of /z-transform, but in terms of duality between two families of orthogonal polynomials in dimension 2. It shows 
that the law of a Brownian motion in the plane, conditioned not to reach the boundaries of an equilateral triangle, 
has the law of the spectrum of an Brownian SU (3) matrix. 

This example, closely related to root systems and reflection groups in the plane, consists in observing the image 
of a planar Brownian motion reflected along the edges of an equilateral triangle. This triangle generates a triangular 
lattice in the plane, and this image is observed through some function Z: R 2 i—»• R 2 which has the property that any 
function R 2 i—»• R which is invariant under the symmetries among the lines of the lattice is a function of Z. This 
image of R 2 through the function Z is a bounded domain in R 2 , with boundary the Steiner’s hypocycloid. 

The Steiner hypocycloid (also called deltoid curve) is the curve obtained in the plane by rotating (from inside) a 
circle with radius 1 on a circle with radius 3. Is is the boundary of a bounded open region in the plane which we call 
the deltoid domain Qq. It is an algebraic curve of degree 4. It’s equation may be written in complex coordinates as 
{D(Z,Z) = 0}, where D is defined in Proposition 5. 



Fig. 1.1 The deltoid domain. 

Consider the following application R 2 i—>• R 2 , which is defined as follows. Let ( 1 , 7 , 7 ) be the three third roots 
of units in the complex plane C, and, identifying R 2 with C, let Z(z ): R 2 1—1 R 2 be the function 

Z(z) = ^(exp(i(l -z)) + exp(j( 7 . ■z)) + exp(i(j-z))'j, 

where z\ ■ Zi denotes the scalar product in R 2 . 

We have 

Proposition3. Let L be the lattice generated in the plane by the points with coordinates M\ = (0,4tt/3) and 
M 2 = (2?r/3,2 n/s/3), and T the (equilateral) triangle with edges {(0,0),Mi,M 2 }. 
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1. The image o/R 2 under the function Z is the closure Qo of the deltoid domain. 

2. Z : R 2 K > R 2 is invariant under the symmetries along all the lines of the lattice L. Moreover, it is injective on the 
triangle T. 

We shall not give a proof of this, which may be checked directly. We refer to [32] for details. As a consequence, 
any measurable function R 2 i—>• R which is invariant under the symmetries of L may be written /(Z), for some 
measurable function / : L2d H > R. 

The particular choice of this function Z is due to the fact that the Laplace operator in R 2 has a nice expression 
through it. Using complex coordinates as described in Section 1.2, one has 

Proposition 4. For the Laplace operator A in R 2 and its associated square field operator T, one has 

rr(z,z)=z-z 2 ,r(z,z) = z-z 2 , 

< F(Z,Z) = 1/2(1-ZZ), (1.8) 

I A(Z) = -Z,A(Z) = -Z, 

This may be checked directly. One sees that the Laplace operator in R 2 has an image through Z in the sense 
described in Section 1.2, given in Proposition 4. This describes the generator of the Brownian motion in the plane, 
reflected along the lines of this lattice, coded through this change of variables. One may express the image measure 
of the Lebesgue measure on the triangle in this system of coordinates. With the help of formula (1.5), we get 

Proposition 5. Let D(Z,Z) = F(Z,Z) 2 — F(Z,Z)F (Z,Z), where F is given by equation (1.8). Then, 

1. D(Z,Z) is positive on Q/>. 

2. {Z3(Z,Z) = 0} is the deltoid curve (that is the boundary of LI d). 

3. The reversible measure for the image operator described by (1.8) has density Z)(Z,Z)~ 1,/2 with respect to the 
Lebesgue measure. 

4. If we write z i = exp(/(l -z)), zi = exp Z3 = exp(i'(y -z)), then 

D(Z,Z ) = -(zi-Z2) 2 (Z2-Z 3 ) 2 (Z3-Zl) 2 /(2 2 3 3 ). 

Remark 3. Observe that thanks to the fact that |z,| = 1 and Z1Z2Z3 = L the expression (zi — Z2) 2 (z2 — Z 3 ) 2 (Z3 — zi) 2 
is always non positive. Moreover, given a complex number Z in the deltoid domain LId, there exist three different 
complex numbers (zi,Z2,Z3) with |z,-| = 1 and Z1Z2Z3 = 1 such that Z = |(zi +Z2 + z 3 )- They are unique up to 
permutation, and are the solutions of X 3 — 3ZX 2 + 3ZX —1=0. Indeed, for such numbers zi,Z2,Z3, 

1 1 1 

3Z — Z1+Z2FZ3 —-1-1-— Z2Z3 + Z1Z3 + ziZ2- 

Zl Z2 Z3 

One may now consider the family of operator * defined through 

|F(Z,Z)=Z-Z 2 ,F(Z,Z)=Z —z 2 , 

lr(Z,Z) = 1/2(1 -ZZ), (1.9) 

[SfW(Z) = -XZ,J?W(Z) = -A Z, 

which is symmetric with respect to the measure = D(Z,Z) t ' 2 ^^/ 6 dZ, with support the set {D(Z,Z) > 0} 
(where dZ is a short hand for the Lebesgue measure in the complex plane) as a direct (although a bit tedious) 



computation shows from a direct application of formula (1.5) (see Section 1.4 for a proof in a general context 
which applies in particular here). 

This family of operators plays a role similar in this context to the one played by the family jP a ,\3 for Jacobi 
polynomials introduced in Section 1.3.2 or for the family S?( a ) introduced in Section 1.3.3 for Laguerre polyno¬ 
mials. 

This density equation (1.5) indicates that, for any pair of smooth functions compactly supported in { I)(Z.Z) > 
0}, the integration by parts (1.4) holds true. Indeed, we have a much stronger result, which extends this formula 
to any pair of smooth functions defined in a neighborhood of Q. This relies of some miraculous property of r)Q 
itself, which has as boundary equation {D(Z,Z) = 0} and for which 

{r{Z,Z)d z D + r{Z,Z)d 2 D = ~3ZD, 

[r(z,z)d z D+r{z,z)d 2 D = ~3ZD. 11 

In particular, r(Z,D) and r(Z,D) vanish on { I) = 0}. This is a sufficient (and indeed necessary) for the integration 
by parts formula (1.4) to be valid for any pair smooth functions restricted on the set {D > 0}, in particular for 
any pair of polynomials (see [4]). Since on the other hand the operator maps polynomials in (Z,Z) into 
polynomials, without increasing their total degrees, the restriction of 2z?^ on the finite dimensional space of 
polynomials with total degree less than k is a symmetric operator (with respect to the L 2 (/^-Euclidean structure) 
on this linear space. We may therefore find an orthonormal basis of such polynomials which are eigenvectors for 
and therefore construct a full orthonormal basis of polynomials made of eigenvectors for S£ { ’ L) . 

These polynomials are an example of Jack’s polynomials associated with root systems (here the root system Az), 
see [24, 15], generalized by MacDonald [23, 24, 25], see also [17, 18, 8], and for which the associated generators 
are Dunkl operators of various kinds, see [21, 20, 14, 28, 30]. 

For A = 4, it turns out that this operator is, up to a scaling factor 8/3, the image of the Laplace (Casimir) 
operator on SU (3) acting on the trace of the matrix. More precisely, on the compact semi-simple Lie group SU (3), 
we associate to each element E in the Lie algebra Sf a (right) vector field A), as follows 

X E (f)(g) = d t (f( g e tE )) lt=0 . 

Then, one choses in the Lie algebra Sf an orthonormal basis E, for the Killing form (which is negative definite), and 
we consider the operator S£ = This is the canonical Laplace operator on the Lie group, and it commutes 

with the group action, from left and right : if L g (f)(x) = f(xg ), and R g (f)(x) = f(gx), then S?L g = LgSS and 
S?R g = R g SS. For the Casimir operator acting on the entries (zij ) of an SU (cl) matrix, one may compute explicitly 
this operator, and obtain, up to a factor 2, the following formulae 

r S£ su ^\z kl ) = — 2 {d ~ l Y +l) z k ,, S? su W{z k i) = 

< (1.11) 
^ ( ! (l (z,kl , Zrq ) — 2 Z k qZrl T fjZk/Zrq - L (z k l , Z n j) = 2.(8 kr 8/q rfZ k lZrq ) * 


A Brownian motion on SU (d) is a diffusion process which has this Casimir operator as generator (there are of 
course many other equivalent definitions of this Brownian motion). 

On SU(3), if one considers the function SU(3) i-»- C which to g € SU(3) associates Z(g) = ^trace(g), then 
one gets for this function Z and for this Casimir operator, an image operator which is the operator |Sf( 4 \ where 
SfW is defined through equation (1.9). Of course, one may perform the computation directly, or use the method 
described in paragraph 1.5.2 to compute from the operator given of SU(d) through formulas (1.11), the actions 
of the generator and the carre du champ on the characteristic polynomial P(X) = dctfATd g) (see also [6] for 
another approach, together with [7] for nice connections with the Riemann-Zeta function). 
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It is worth to observe that functions on SU (3) which depend only on this renormalized trace Z are nothing else 
but spectral functions. Indeed, if a matrix g £ SU{3) have eigenvalues (Ai, A-2,-^-3), with |A,j = 1 and A1A2A3 = 1 , 
then a spectral function, that is a symmetric function of (Ai. A2, A3), depends only on Aj +A2 + A3 = 3Z and, as 
observed in Remark 3, AiA2 + A2A3 + A3A1 = 3Z. 

Then, using the function D which is the determinant of the metric involved in equation f 1 . 10 ), one may check 
directly that 

_S?W(£>(Z,Z) (5-2;i)/6 ) = (2A—5 )£>(Z,Z) (5 ~ 2A)/6 , 

so that one may use the function h = DiZ.Zf 5 24 R 6 to perform an h transform on and we obtain 


Indeed, as we shall see in Section 1.4, this /z-transform identity relies only on equation (1.10). In particular, 
moving back to the triangle through the inverse function Z 1 , for A = 1, which corresponds to the Brownian 
motion reflected at the boundaries of the triangular lattice, the h transform is ..Z 441 , which corresponds to the 
spectral measure on SU (3). Then, for this particular case A = 1, we get 

Proposition 6. A Brownian motion in the equilateral triangle T, conditioned to never reach the boundary of the 
triangle, has the law of the image under Z” 1 of the spectrum of an SU (3) Brownian matrix 


1.3.5 An example in the unit ball in W 1 

Another example comes from the spherical Brownian motion on the unit sphere 

§ d = {(xi,---,x d+ 1 )eR d+ 1 ,£x 2 = i}. 

i 

To describe the Brownian motion on S d , we look at its generator, that is this the spherical Laplace operator may. It 
may be described through its action on the restriction to the sphere of the coordinates x h seen as functions E> d i-»- R. 
Then.for the Laplace operator A “ and its associated carre du champ operator F, one has 

(xi) = —dxi, r(xj,Xj) = Sjj—XiXj. ( 1 . 12 ) 

This operator is invariant under the rotations of M 4 11 , and as a consequence its reversible probability measure 
is the uniform measure on the sphere (normalized to be a probability). A system of coordinates for the upper 
half sphere {xd+ 1 > 0} is given by (xi,--» ,x^) £ By, where B^ = jXfxf = ||x|| 2 < 1} is the unit ball in R rf . In 
this system of coordinates, and thanks to formula (1.5), one checks easily that, up to a normalizing constant, the 
reversible measure is (1 — ||x|| 2 )~ l / 2 <ix, which is therefore the density of the uniform measure on the sphere in this 
system of coordinates (see [2]). 

Now, one may consider some larger dimension m > d and project the Brownian motion on on the unit ball in 
R (/ through (xi, • • • ,x m +i) 1 —> (xi, • • • ,xf). Formula (1.12) provides immediately that this image is again a diffusion 
process with generator 

J2? (M) (x/) = -mxj,r(xi,xj) = Sij-XjXj, (L13) 

that is the same formula as (1.12) except that now m is no longer the dimension of the ball. Once again, 
formula (1.5) provides the reversible measure for this operator, which is, up to a normalizing constant, (1 — 
||x|| 2 j(m-i d )/ 2 dx, which is therefore the image measure of the uniform measure of the sphere through this projec¬ 
tion. 



As before, the boundary of the domain (the unit ball) has equation {1 — ||x|| 2 = 0}, and we have a boundary 
equation 

r(x h log(l - ||x|| 2 )) =-2x«, (1.14) 

similar to equation (1.10). 

Now, it is again easily checked that, for the function h = (1 — \\x\\ 2 )~^ m ~ l ~ d ^ 2 , one has 


^ m \h)=d{m-d- 1 )h, 


so that one may perform the associated /z-transform for which 


(^>(m))(fc) _ j£>(2d+2- 


In the case where m = d, on sees that , ( Z Ad] , which is the Laplace operator in this system of coordinates, is trans¬ 
formed into _S ?( d + 2 ), which is the projection of the spherical Laplace operator in Z 1 ' 2 onto the unit ball in 1RJ? 

As a consequence, we get 


Proposition 7. A spherical Brownian motion on the unit sphere S d C R </+1 conditioned to remain in a half sphere 
{xd +1 > 0}, has the law of the projection of a spherical Brownian motion on S d+2 onto the unit ball in R rf , lifted 
on the half upper sphere in R^ +1 . 


1.4 General /i-transform for models associated with orthogonal polynomials 


We shall see in this section that all the above examples appear as particular examples, or limit examples, of a 
very generic one when orthogonal polynomials come into play. Everything relies on a boundary equation similar 
to (1.10) or (1.14), which appears as soon as one has a family of orthogonal polynomials which are eigenvectors 
of diffusion operators. 

Let us recall some basic facts about diffusion associated with orthogonal polynomials, following [4], We are 
interested in bounded open sets £2 C R d , with piecewise Z ' 1 boundary. On £2, we have a probability measure 
/i with smooth density p with respect to the Lebesgue measure, and an elliptic diffusion operator Jz? which is 
symmetric in L 2 (p). We suppose moreover that polynomials belong to the domain of Jz?, and that Jz? maps the set 
of polynomials with total degree less than k into itself. Then, we may find a L 2 (p) orthonormal basis formed 
with polynomials which are eigenvectors for Jz?. Following [2], this is entirely described by the triple (£2.T .p), 
where F is the square field operator of Jz?. 

We call such a system (f2,F,p) a polynomial system. 

Then, one of the main results of [4] is the following 

Theorem 1. 

1. The boundary d£2 is included in an algebraic surface with reduced equation {P = 0}, where P is a polynomial 
which may we written as P\ ■ ■ ■ where the polynomials Pi are real, and complex irreducible. 

2. If Jz? = fjijg' 1 df: + Y,i b‘d h where the coefficients g’i are degree at most 2 polynomials and b’ are polynomials 
with degree at most 1. 

3. The polynomial P divides det(g'- / ) (that we write det(F) in what follows, and which is a polynomial with degree 
at most 2d). 

4. For each irreducible polynomial P r appearing in the equation of the boundary, there exist polynomials L\ r with 
degree at most 1 such that 
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Vi=l,«--,d, '£g ij dlogP r = L ir . (1.15) 

j 

5. Let £2 be a bounded set, with boundary described by a reduced polynomial equation {P] ■ ■ ■ !\ = 0}, such that 
there exist a solution (g lJ ,Li k ) to equation (1.15) with (gpositive definite in £2. Call T(f,f ) = Y,ij 8 ' J difdjf 
the associated squared field operator. Then for any choice of real numbers {CC\, ■ ■ ■ , Ct k } such that P ® 1 ■ ■ ■ P k ' k is 
integrable over £2 for the Lebesgue measure, setting 

,■■■ ,a k (dx) = Ca j ,a k P\ 1 • • • P^ k dx, 

where C ai ,-,a k is a normalizing constant, then (£2 , p ai ... tCCk ) is a polynomial system. 

6 . When P = Cdet(F), that is when those 2 polynomials have the same degree, then there are no other measures 
p for which (£ 2 ,r,p) is a polynomial system. 

Remark 4. Equation (1.15), that we shall call the boundary equation (not to be confused with the equation of the 
boundary), may be written in a more compact form F(x,-. log/)-) = L ir . Thanks to the fact that each polynomial P r 
is irreducible, this is also equivalent to the fact that F(x,\ logP) = L,, for a family L, of polynomials with degree at 
most 1. 

One must be a bit careful about the reduced equation of the boundary {P = 0}, when P = P\ ■ ■■ . This means 
that each regular point of the boundary is contained in exactly one of the algebraic surfaces {P-fx) = 0}, and that 
for each / = 1 • • • k, there is at least one regular point x of the boundary such that P, (x) = 0. In particular, for a regular 
point x £ d£2 such that Pj(x) = 0, then for j zfi i, Pj(x) f 0 in a neighborhood of such a point, and P,(x) = 0 
in '2/ D r)£2. It is not too hard to see that such a polynomial If if real irreducible, is also complex irreducible (if 
not, it would be written as P 2 + Q 2 , and P = Q = 0 on "2/ n c)£2). It is worth to observe that since P divides det(F) 
and that (g‘ J ) is positive definite on £2, then no one of the polynomials I) appearing in the boundary equation may 
vanish in £2. We may therefore chose them to be all positive on £2. 

The reader should also be aware that equation (1.15), or more precisely the compact form given in Remark 4, 
and which is the generalization of equations (1.10) and (1.14), is a very strong constraint on the polynomial P. 
Indeed, given P, if one wants to determine the coefficients (g lJ ) and L,, this equation is a linear equation in terms of 
the coefficients of g l l and L\, for which we expect to find some non vanishing solution. But the number of equations 
is much bigger than the number of unknowns, and indeed very few polynomials P may satisfy those constraints. 
In dimension 2 for example, up to affine invariance, there are exactly 10 such polynomials, plus one one parameter 
family (see [4]). The deltoid curve of paragraph 1.3.4 is just one of them. 

Remark 5. We shall not use the full strength of this theorem in the examples developed here. The important fact 
is the boundary equation (1.15), which may be checked directly on many examples, and is the unique property 
required for the general /z-transform described in Theorem 2. 

Given a bounded set £2 and an operator F satisfying the conditions of Theorem 1, and for any choice of 
{ai, • • • , a . k } such that • • • P® k is integrable over £2 for the Lebesgue measure, we have a corresponding sym¬ 
metric operator f£ ax ... a ,. For this operator, as was the case in paragraphs 1.3.4 and 1.3.5, one may extend the 
integration by parts (1.4) to any pair of polynomials, and this provides a sequence of orthogonal polynomials 
which are eigenvectors of the operator L£ ai ,a k - 

Conversely, the boundary equation (1.15) is automatic as soon as we have a generator on a bounded set with 
regular boundary, and a complete system of eigenvectors which are polynomials. But it may happen that those 
conditions are satisfied even on non bounded domains, and even when the associated measure is infinite (this 
appears in general in limits of such polynomial models, as in the Laguerre and Bessel cases). We may therefore 
give a statement in a quite general setting. 



Theorem 2. Assume that a symmetric positive definite matrix (g^) on some open set £2 C W 1 , is such that for any 
( i,j ), g IJ is a polynomial of degree at most 2. Let us call T the associated square field operator. Suppose moreover 
that we have some polynomials P k , positive on £2, such that, for any k, 

Vi = l,"• ,d, Y^g ij djlogPr = £r(x ( ',log/\) = L ik , (1.16) 

i i 

where Lj^ are degree 1 polynomials. For any («i, • - - ,CC k ), let p ai .--,a k be the measure with density P■ ■ ■ P k k 
with respect to the Lebesgue measure on £ 2, and let Jz?<xi,-,a k be the generator associated with the Markov triple 

{£2,r,p a i,a*)- 

Then, there exist constants c k such that, for any (ai, • • • , a k ), the function h = P l ai ■ ■ ■ P k satisfies 

■%a w -,a k {h) = ~(^a k C k )h. 

k 


Moreover, -,a k ) {ll) = ff-a u - ,-a r 

Proof. We shall prove the assertion with c k = £,• diLi.k- 
With p = P® 1 • ■ ■ P k k , we write our operator S£a x ,--,a k as 

Z^i+Z^ 

ij l 


where 

b, = Y djg lJ + Y a rS ,Jd j lo g P r = Y d JS ,j + Y a ' L i.r- ( 1 - 17 ) 

j r,j j r 

With 

ij i 

then 

^o!,-.«*(/) = ^o(/)+£«/r(iogP,-,/). (i.i8) 

i 

What we want to show is (h) = ch, or 

,a k {\ogh) + r(logh,\ogh) = c. 

With, log/; — a, log Pj, and comparing with equation (1.18), this amounts to 

-^o(log h) = -Y a ‘i-%’o(l°gPr) = c- 


We may first take derivative in equation (1.16) with respect to x, and add the results in i to get 

Yg' J di j log Pr + Y di C? ,; ) dj lo S P< = Y d ‘ L hr = C r > 

ij i i 

that is Jz?o(logPr) = c r . 

It remains to add these identities over r to get the required result. 

Comparing the reversible measures, it is then immediate to check that (jSf ai) ... )Qj . )W = J Sf- a u~,-a k > D 
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Remark 6. The function h is always the inverse of the density with respect to the Lebesgue measure, in the system 
of coordinates in which we have this polynomial structure. Of course, the choice of the coordinate system is 
related to the fact that, in those coordinates, we have orthogonal polynomials (at least when the measure is finite 
on a bounded set). In the Bessel case, for example, which is a limit of a Laguerre models, one has to change x to x 2 
to get a simple correspondance between the h function and the density. The same is true in many natural examples, 
where one has to perform some change of variable to get the right representation (for example from the triangle to 
the deltoid in paragraph 1.3.4). 

Remark 7. In many situations, there are natural geometric interpretations for these polynomial models when the 
parameters (0£i,--- ,a k ) are half integers, in general with a, > —1/2. The case a,- = —1/2 often corresponds to 
Laplace operators, while the dual case a, = 1/2 often corresponds to the projection of a Laplace operator in larger 
dimension. 


1.5 Further examples 


We shall provide two more examples, one which follows directly from Theorem 2, and another one on a non 
bounded domain with infinite measure. One may provide a lot of such examples, many of them arising from Lie 
group theory, Dunkl operators, random matrices, etc. However, we chose to present those two cases because they 
put forward some specific features of diffusion operators associated with orthogonal polynomials. 


1.5.1 Matrix Jacobi processes 


This model had been introduced by Y. Doumerc in his thesis [13], and had also been studied in the complex case, 
especially from the asymptotic point of view in [9, 10]. It plays a similar role than the one-dimensional Jacobi 
processes for matrices. One starts from the Brownian motion on the group SO(d). Since SO(d) is a semi-simple 
compact Lie group, it has a canonical Casimir operator similar to the one described in equation (1.11). If O = (my) 
is an SO(d) matrix, then the Casimir operator may be described through it’s action on the entries my. One gets 

if(my) = -(c/- l)my, r(m k i,m qp ) = 8^ qp ) -m kp m q i. (1.19) 

Observe that when restricted to a single line or column, one recovers the spherical Laplace operator on Ef l 1 
described in equation (1.12). 

An .S'0(c/)-Brownian matrix is then a diffusion process with generator this Casimir operator on SO{d). 

It is again clear from the form of the operator ££ that it preserves for each k G N the set of polynomials in the 
entries (my) with total degree less that k. However, these "coordinates" (my) are not independent, since they satisfy 
algebraic relations, encoded in the fact that 00* = Id. We may not apply directly our main result Theorem 2. We 
shall nevertheless look at some projected models on which the method applies. 

One may extract some pxq submatrix N by selecting p lines and q columns, and we observe that the generator 
acting on the entries of this extracted matrix N depend only on the entries of N. Therefore, the operator projects on 
these extracted p x q matrices and the associated process is again a diffusion process : we call this the projection 
of the Brownian motion in SO(d) onto the set ./M vq of p x q matrices. Thanks to formula (1.5), one may compute 
the density of the image measure, with respect to the Lebesgue measure in the entries of N. Whenever p + q <d, 
it happens to be, up to a normalizing constant det(Id — NN*)( d ~ l ~ p ~ q ^ 2 , with support the set Q = {N,NN* < Id}. 
This formula is easy to check if we recall that, for a matrix M with entries (my), 



log det(M) = Mjf 1 , 


a consequence of Cramer’s formula. 

When p + q > d + 1, there are however algebraic relations between the entries of N and the image measure has 
no density with respect to the Lebesgue measure. For example, when p + q = d + 1, then the measure concentrates 
on the algebraic set {det(ld — AW*) = 0}. It may be checked that it has a density with respect of the Lebesgue 
measure of this hypersurface. Indeed, one may fix p and q and consider d as a parameter. It is worth to observe that 
the function det(Id — NN*) a is not integrable on the domain £1 whenever a < — 1. Moreover,, when a > — 1 and 
a —> — 1, the probability measure with density C a det(Id — NN*) a concentrates on the set {det(Id — AW*) = 0}, 
and the limit is a measure supported by this surface with a density with respect of the surface measure. Things 
become even worse as the number p + q increases, the measure being concentrated on manifolds with higher and 
higher co-dimensions. 

We are in a situation different from the sphere case here, since we may not chose the parameters in which 
the operator has a nice polynomial expression as a local system of coordinates. Indeed, the Lie group SO(d) is 
a d{d — l)/2 manifold. Since we want algebraically independent coordinates, we are limited to pq ones, with 
p + q < d, we may have at most d 2 /4 algebraically independent such polynomial coordinates, which for d > 2 is 
less than the dimension of the manifold. 

It is worth to observe that, again when p + q < d, one has pq variables, the determinant of the metric F is a 
degree 2pq polynomial, whereas detfld - NN*) = detfld — N*N) is of degree at most 2 min (p,q). We are not in the 
case of maximal degree for the boundary equation. When p + q — d, the density measure is detfld-AW*) ' ,2 , but 
the corresponding operator is not a Laplace operator (for which the density of the measure would be det(F) -1 / 2 ). 
Since we are in the situation of orthogonal polynomials as described in Section 1.4, we know that we may perform 
an /z-transform. 

For the particular case where d = p + q, we get 

Proposition 8. The matrix N projected from an SO(d)-Brownian matrix on conditioned to remain in the set 
{AW* < Id} has the law of the projection of a SO(d + 2)-Brownian matrix on 


1.5.2 Brownian motion in a Weyl chamber 

This last example is again quite well known, but is happens to fit also with the general picture associated with 
orthogonal polynomials, although no orthogonal polynomials are associated with it. Indeed, it does not follow 
directly from the setting of Section 1.4, one the one side because it is non compact, on the other because the 
reversible measure in this situation is infinite. But it satisfies the all the algebraic properties described is Section 1.4, 
and we may then check that we may apply the result for the associated /z-transforms. Indeed, one may replace in 
what follows Brownian motion by Ornstein-Uhlenbeck operators, which have as reversible measure a Gaussian 
measure with variance cr 2 , and then let a go to infinity. In the Ornstein-Uhlenbeck case, we are in the setting 
of orthogonal polynomials, however with a non bounded domain. But this would introduce further complication, 
since the Brownian case gives simpler formulas. 

As described above, the /z-transform is easy to compute in a system of coordinates which have some relevant 
polynomial structure. Here, one good choice for the coordinate system are the elementary symmetric functions in 
d variables. We shall perform mainly computations on these elementary symmetric functions of the components of 
the d -dimensional Brownian motion, following [5], In 9 d , one may consider the Brownian ( B] ,■■■ .Bf) and reflect 
it around the hyperplanes which are the boundaries of the set {x\ < ■■■ < x,j }, which is usually called a Weyl 
chamber. To describe this reflected Brownian motion, it is easier to consider the elementary symmetric functions 
which are the coefficients of the polynomial 
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P(X) = f\(X~x i ) = Y j a i X i , 
i= 1 i=0 

where a d = 1 and the functions aj,i = 0 , ■ • -d — 1 are, up to a sign, the elementary symmetric functions of the 
variables (x,). The map (x,j i—>• (a,-) is a diffeomorphism in the Weyl chamber {xi < • • • < x^} onto it’s image. To 
understand the image, one has to consider the discriminant disc(f’), a polynomial in the variables (a,), which is, 
up to a sign (—1 ) d ( d ~P/ 2 , the following ( 2 d — 1) x ( 2 d — 1) determinant 
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2.02 a 1 / 


It turns out that this discriminant is Yh<j( x j ~ x i) 2 - The image of the Weyl chamber is the connected component 
£2 of the set {disc(f’) 0 } which contains the image of any polynomial with d real distinct roots, and the image of 
the boundary of the Weyl chamber is d£2, a subset of the algebraic surface {disc(f’) = 0}. It is not hard to observe 
(by induction on the dimension d) that the image of the Lebesgue measure dx\ ■ ■ ■ dx t j on the Weyl chamber is 
nothing else than lfldisc(/ > )“ I / 2 n^ fl i- 

Now, the Brownian motion in R' y may be described, up to a factor 2, through 

r(x h Xj) = 8jj, A(xj) = 0. 

We want to describe this operator acting on the variables (a o, • • • , a ( \ i). Since any of the functions aj is a degree 
1 polynomial in the variables x/, one has A (aj) = 0, j = 0, • • • , d. To compute F(a,, ay), it is simpler to compute 

r(P(X),P(Y))=Y,X i Y l r(a i , aj ). 

hj 


We obtain 

Proposition 9. The image of the operator A in R" on the coefficients of the polynomial P(X) = — xf) is given 

by 

r(P(X),P(Y)) = v L l (P'(X)P(Y)-P'(Y)P(X)), A(P(X))= 0. (1.20) 

Proof The second formula is a direct consequence of A(af) = 0, while for the first, it is simpler to look at 
r(logP(X),logP(F)). 

r(iogP(x),io g p(y)) = £r(iog(z-x«),iog(F-x;)) 
ij 

= L (X _J ( y_ Xj ) r (^>) 

_ y i _ i (pyn pyp\ 

j (X — xf)(Y — xf) Y-X\P(X) P(Y) ) 










□ 


Remark 8. From formula (1.20), it is clear that r{a-,.aj) are degree 2 polynomials in the variables a,-. 

The image of the Brownian motion B t in the variables (a,-) is nothing else than the Brownian motion reflected 
through the walls of the Weyl chamber. Its generator is described through the F operator given in equation (1.20) 
and it is the image of the Laplace operator on the Weyl chamber. Since it is an Euclidean Laplace operator, the 
reversible measure is, up to a constant, det(F) l 2 , and this shows that the determinant det(F) of the metric is, up 
to a constant, disc(P). 

Moreover, from the general representation of diffusion operators (1.3), and the equation (1.5) giving the re¬ 
versible measure, we have, with p = disc(F)" 1 / 2 , bj = 0, 

Y / r(a i ,a j )d ai \ogp = (1.21) 

ij j 

Since d aj r{ai 1 af) is a degree at most one polynomial in the variables a;, this is nothing else than the boundary 
equation (1.15) for general polynomial models. We may therefore apply the general result described in Section 1.4. 

In order to identify the result of the /z-transform, an important formula relating F and the discriminant function 
is the following 

Proposition 10. For the operator F defined in (1.20), one gas 

F(P(X),logdisc(P)) = -P"(X). (1.22) 


Proof. One may find a proof of this formula in [5], but the one we propose here is simpler. To check equation (1.22), 
it is enough to establish it it in a Weyl chamber (xi < x\ < ■ ■ ■ < xj} where P(X) = n(X ~ x i) and disc(P) = 
n,<,(x, — X j)2, since the map (x\ ,• • • ,x^) H > P(X) is a local diffeormorphismin this domain. 

In those coordinates, F (xj.Xj) = 8, ; and, from the change of variable formula (1.2), one has 


F (log P (X), disc (P)) =2 Y, F(log (X — x,), log (xj —x^)) = -2 £ — - -—(Sjj - 8 ik ). 

i,j<k i,j<k A X ' X J Xk 

From which one gets 

F(logP(X),logdisc(P)) = -2£-i-—. 

faX-XiXi-Xj 

On the other hand. 


P" 

~P 


/P'y (P '\2 „ l 

\~PJ + \~PJ -^.(X-x^X-xj) 

1 V y , V y , ^ y , y . ^ ~V^ y , y , y , 

• / j A A-i A Ay A; A j j S\- A; A; A j 


From this we get 


which in turns gives (1.22). □ 


F(logP,logdisc(P)) = -y, 


Proposition 10 is central in the identification of various processes with the same F given by (1.20). It turns out 
that the same operator with this F operator and reversible measure disc(P) 1/ 2 has a nice geometric interpretation: 
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namely, it is the Dyson complex process, that is the law of the spectrum of Hermitian Brownian matrices, intro¬ 
duced by Dyson [16]. In the same way, the case where the reversible measure is the Lebesgue measure corresponds 
to Dyson process for real symmetric matrices, and p = disc(B) 3 / 2 corresponds to Dyson process for symmetric 
quaternionic matrices , see [1, 5, 19]. 

Let us show a direct way to check this (first in the real symmetric case, where it is simpler). The Brownian 
motion on symmetric matrices is nothing else that the Brownian motion of the Euclidean space of symmetric 
matrices M , endowed with the Euclidean norm ||M|| 2 = trace (M 2 ). When M = (m,y), this may be described as 


r(mij,m k i) = ^{8 ik 8ji + 8u8 jk ), -S?(m /; -) = 0. 

One may look at its action of the characteristic polynomial P(X) = detfATd — M). We get 
Proposition 11. For the characteristic polynomial associated with a Brownian symmetric matrix, one has 

rqcg/TO.iog/W) = y±x(m-m)’ = 

Proof. To compute r(P(X),P(Y)) and 2 z?(P(A)). In order to apply the change of variable formula (1.3), we may 
apply the general formulas for the determinant function 

d„ Hj logdet M = MJ j \ d mjj d,„ kl logdet M = -Mf k M n 1 , 

which are direct consequences of Cramer’s formulas for the inverse matrix. 

Then the formulas are direct applications of the chain rule formula. □ 

We may now compare with equation (1.22) to see that the reversible measure for the spectral measure for 
Brownian symmetric matrices, given by the general formula (1.5), in the system of coordinates which are the 
coefficients (a,) of the characteristic polynomial, is the Lebesgue measure. 

We may perform the same computation for Hermitian matrices. In this situation, one would consider a complex 
valued matrix M with entries (zi j ) and satisfying 

r(zij,zici) = o, r(zij,Zki) = 8 n 8j k , ££{zij) = 0. 

One may again perform the same computation on P(X) = det(XId — M), and we get 
Proposition 12. For the characteristic polynomial associated with a Brownian Hermitian matrix, one has 

F(P(X),P(V)) = (P'(X)P(Y) - F f (Y)P(X)) , F£P(X) = -P". 

We do not give the proof, which follows along the same lines that the one of Proposition 11. More details may 
be found in [5], 

As a consequence, comparing with equation (1.22) and equation (1.5) in the system of coordinates given by the 
coefficients of P(X), the density of the reversible measure for the Hermitian Dyson process is disc]/ J ) 1/ 2 whereas 
the density of the reversible measure of the Brownian motion in the Weyl chamber is disc(P) -1//2 . 

Transfering back to the Weyl Chamber through the local diffeomorphism between the coefficients of P(X) and 
the roots (xi < xo < ■ ■ ■ < xf) of P(X). We obtain 

Proposition 13. The Brownian motion conditioned not to reach the boundary of the Weyl chamber {v| < • • • < xj} 
has the law of the spectrum of an Hermitian d x d matrix. 
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